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Abstract 

We consider the current experimental constraints on the parameter space of the MSSM and 
NMSSM. Then in the allowed parameter space we examine the Higgs pair production at the 14 
TeV LHC via bb > hh (h is the 125 GeV SM-like Higg boson) with one-loop SUSY QCD correction 
and compare it with the production via gg — hh. We obtain the following observations: (i) For 
the MSSM the production rate of bb — hh can reach 50 fb and thus can be competitive with 
gg — hh, while for the NMSSM bb — hh has a much smaller rate than gg — hh due to the 
suppression of the hbb coupling; (ii) The SUSY-QCD correction to bb — hh is sizable, which can 
reach 45% for the MSSM and 15% for the NMSSM within the lo region of the Higgs data; (iii) In 
the heavy SUSY limit (all soft mass parameters become heavy), the SUSY effects decouple rather 
slowly from the Higgs pair production (especially the gg — hh process), which, for Msusy = 5 
TeV and ma < 1 TeV, can enhance the production rate by a factor of 1.5 and 1.3 for the MSSM 
and NMSSM, respectively. So, the Higgs pair production may be helpful for unraveling the effects 


of heavy SUSY. 


PACS numbers: 


I. INTRODUCTION 


The discovery of a Higgs boson at around 125 GeV has been announced by the ATLAS 
and CMS collaborations il Up to now, the measurements of the Higgs boson properties 
are in good agreement with the Standard Model (SM) predictions except for the enhanced 
diphoton rate c/o sy = 1.651636 reported by the ATLAS collaboration. The future precise 
measurements will further test the SM and allow for a probe for new physics like super- 
symmetry (SUSY) which is a promising framework to accommodate such a 125 GeV Higgs 
boson |a. Therefore, the intensive studies of the Higgs productions and decays are very 
important and urgent. 

Among the productions of the Higgs boson at the LHC, the pair production is a rare 
process but quite important since it can be used to measure the Higgs self-couplings 
On the experimental side, the discovery potential of Higgs pair signal at the LHC has 
been studied by analyzing the decay channels hh > bbyy/bbu* u- q Recently, the li 


substructure technique was applied to the Higgs pair production in the boosted final states 


p] 


such as hh — bbr*r- /bbW +W , which was found to be powerful in observing the 
events at the 14 TeV LHC with 600 fb~! integrated luminosity ld. On the theoretical side, 
in the SM the main pair production mechanism is found to be the gluon fusion gg — hh 
via heavy quark loops . Numerous studies have also been performed for Higgs pair 
production in new physics models . Note that although the bottom quark annilation 

bi bÀ. it can be 
significantly enhanced via the enlarged Abb coupling in new physics models like the Minimal 


Supersymmetric Standard Model (MSSM) [25]. 


bb — hh has a much smaller rate than the gluon fusion process in the SM 


In this work, we revisit the Higgs pair production in SUSY for two reasons. One is that the 
sizable SUSY-QCD correction must be considered for bb — hh, which has been presented in 
the MSSM but not in the NMSSM bd bj. The other is that the studies should be updated 
by using the latest experimental constraints including the recent LHC Higgs data, the LHCb 
B, — wtp” data and the Planck dark matter relic density. It is also notable that the masses 
of the third generation sparticles involved in the SUSY-QCD correction to bb — hh have 
been pushed up to a few hundred GeV by the LHC direct searches . So the size of such 
a correction will be quite different from the previous results in the literature 25) f 


This paper is organized as follows. In Section II we briefly review the Higgs sectors in the 


MSSM and NMSSM and give a description of the analytic calculation of the SUSY-QCD 
correction. Then in Section III we present the numerical results of Higgs pair production at 
the LHC and discuss the SUSY-QCD residual effects in the heavy sparticle limit. Finally, 


we draw the conclusion in Section IV. 


II. A DESCRIPTION OF MODELS AND ANALYTIC CALCULATIONS 


In the MSSM there are two complex Higgs doublets, H, and Hz, which give rise to five 
physical Higgs bosons: two CP-even (h, H), one CP-odd (A) and a charged pair (H*). 


Due to the u term appearing in the superpotential, the MSSM suffers from the u-problem. 
Besides, in order to give a 125 GeV SM-like Higgs boson, large corrections to the Higgs 
mass from heavy stops is needed, which will lead to the little fine tuning problem . To 
overcome these difficulties, we can go beyond the MSSM. One alternative is the NMSSM, 
which introduces a singlet Higgs field. In the NMSSM the u term does not appear in the 
superpotential. Instead, it is generated when the singlet Higgs field develops a vev. Also, 
the SM-like Higgs boson gets an extra tree-level mass from the mixing with the singlet field 
and thus the stops are not necessarily heavy to push up the Higgs mass, which alleviates the 
little fine-tuning problem [ad l3. In the NMSSM the singlet Higgs field mixes with the other 
two doublet scalars. Then the Higgs sector contains seven Higgs bosons, i.e., compared with 
the five Higgs bosons in the MSSM, the NMSSM contains one more CP-even and one more 
CP-odd Higgs bosons. In the following H; denote the real scalar components of H,,, in the 
MSSM and Hi denote the real scalar components of H4,,, in the NMSSM. tan f = v/va 
is also used in our paper (here Ha, H,, and H, are the down-type, up-type and singlet Higgs 


fields, respectively). One can get the mass eigenstates from the CP-even states: 


where U2 --U2, = 1, Vi -V24-V2 = 1 and the h; is aligned by mass. The singlet contribution 
is reflected by the rotation matrix elements V;; via the formula hsm = Va ,,1Hi-- Va 4,5 H2 + 
Vi, 3H3 (a large Vi;,,3 means that hsm has a considerable singlet component). 

In our calculations, we follow the simplifed ACOT prescription to deal with the b-quark 
mass . By including the QCD and SUSY-QCD effects to the bottom Yukawa cou- 


plings, we can respectively obtain the effective h;bb couplings in the MSSM |3 and 
NMSSM [46]: 
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where 


y 1 i " 
AMSSM — (rea mr) (i = 1,2), 


1+ A} ? TI tan B 
ANMSSM _ DA h a NE m + a) | (c= 125); 
Ay = = mgzu tan BI (m5, , M$,» ma) 
A e S mgA] (mj L3 m?) Ay = fu (5) 


Here it should be noted that due to the contribution of the singlet field to the effective 


potential, an additional correction term A; a appears in the NMSSM. The vg and v, are 


the VEVs of the Higgs fields H, and Hg respectively. The auxiliary function J is defined as 
1 
(a — b)(b — c)(c — a) 


. The value of mPR is related to the QCD-MS mass mMS (which is usually taken as an input 


I(a,b,c) 2 — (abn $ + bem = + can Ê) (6) 
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parameter ) by 


= (73 — 3n) | , (7) 
where n is the number of active quark flavors and ml? S( ur) is taken as 


— Us(ug, m,)Us(m,,my)my(ms) for pR>m 
m (un) = (8) 


Us (pig, Tip) Tria (ris) for HR <í M. 


When Qə > Q4, the evolution factor U„ reads 


dn — 
Un(Q2, Q1) = (27) h + Oty) esa) ; (9) 
where 
12 8982 — 504n + 40n? 
"= 33-3. 77 7 388-28B C s 


Since the A;-related corrections have already been included into the tree-level contribu- 
tion, we need the following counter terms to subtract them to avoid double counting in the 


one-loop calculations |3 
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For SUSY-QCD corrections to bb — hh, the sbottoms and gluino are involved in the 


loops. The sbottom mass matrix is given by |48 


me ng] 


M? = br : f (13) 
Mp Xp ae 
where 
m2 = m +m- m? E — I aui w) cos(28) 
bo Q ^ 2. 3 ud 
1 
mi. = m 4m- zmz sin? Oy cos(28) , 
Xy = Ay— utan B. (14) 


After diagonalizing Eq.(13), we can obtain the sbottom masses mz, , and the mixing angle 
Oz: 


1 2 
mi. J mi, + mp. F (mj, = m) + 4m2X ] f 
2my X, 
Mi, T Mir 


The Feynman diagrams for one-loop SUSY-QCD corrections to bb + hh has been rep- 
resented in |29]. To preserve supersymmetry, we adopt the dimension reduction method 
to regulate the UV divergences in the gluino and squark loops. Then we use the on-shell 


renormalization scheme to remove these UV divergences. 


IIl. NUMERICAL STUDIES 
A. A scan of parameter space 


We use NMSSM Tools and LoopTools to perform a random scan over the pa- 


rameter space and loop calculations. For simplicity, we assume an universal parameter Mr 


5 


for the slepton sector and fix all irrelevant soft parameters for first two generation of the 
squark sector to be 1 TeV. We also set Mp3 = My3 and A, = A; for the third generation 
of the squarks. Besides, we impose the grand unification relation of the gaugino masses, 
3M, /5a, = Mz/a2 = M3/a3, and treat M; as an input parameter. The parameter ranges 


in our scan are: 
(a) For the MSSM 


1<tan@ < 60, 100 GeV € Ma <1 TeV, 100 GeV < u < 2 TeV 
100 GeV < Mos, Mus € 2 TeV, 100 GeV € Mz3 < 1 TeV 
|A| <5 TeV, 50 GeV € Mı < 500 GeV. (16) 


(b) For the NMSSM 


05€Ax07, 01€ & € 0.51, |A,| € 1 TeV (15) 
1 € tan < 10, 100 GeV < u < 600 GeV, 100 GeV < M4 < 1 TeV 

100 GeV < Mos, Mus € 2 TeV,100 GeV € Mz3 € 1 TeV 

|A| <5 TeV, 50 GeV < Mı < 500 GeV. 


In our scan we consider the following experimental constraints: 


(i) The bounds for Higgs boson from the LEP, Tevatron and LHC experiments and re- 
quire the SM-like Higgs mass to be in the range of 123 GeV « mp < 127 GeV; Here 
we require the surviving samples to explain the observable at 2e level which has an 
experimental central value. For the LEP and Tevatron limits, the upper or lower 


bounds are implemented in our scan. For the LHC Higgs search of H/A — rr|5]] and 


H= > rw,|52J, we require the samples to satisfy the upper limits. 
(ii) The constraints from the precision electroweak data and flavor physics at 2c level; 


(iii) The dark matter relic density from Plank at 3c level and the limit of direct detection 


from XENON100 ; 


(iv) The explanation of muon g — 2 at 2c level [55]. 
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FIG. 1: The measured signal strength of Higgs boson with their lo error-bars and selection effi- 


ciencies €, for each production mode p and decay mode at the 7--8 TeV LHC and Tevatron. 


In our scan, for each experimental data which has a central value, we require the samples 
to agree with the experimental data at 2c level, except for the dark matter relic density 
which is required to agree with the measured value at 30 level (we made such a choice just 
in order to be consistent with the analysis in the literature). For the LEP and Tevatron 


direct search bounds on sparticle masses, we just require the samples to satisfy such bounds. 


For the LHC Higgs search of H/A — rr and H^ — rTv,, we require the samples to satisfy 
the upper limits. The scan ranges of the parameters are large, we keep the samples survived 
various experimental constraints as stated above. Besides, we further require gluino mass 
larger than 1 TeV to avoid multi-jets search on SUSY|56]. However, we did not impose other 
LHC direct limits on sparticles for the following reasons. First, we required the first and 
second generations of squarks to be 1 TeV and the gluino beyond 1 TeV. But the latest LHC 
search results gave more stringent constraints on such squark and gluino mass (the most 
stringent bound is for the CMSSM, which is mg > 1.7 TeV in case of mg ~ m; and mg > 1.1 
TeV in case of mg > mz). Actually, our results are not sensitive to these masses. Second, 
the current LHC limit is about 500-600 GeV for stop and 400-600 GeV for sbottom|57}. 
However, such limits were obtained in some simplified model or by assuming a certain decay 


branching ratio to be 100%. In our case the stop and sbottom decays are quite complicated, 


which will weaken the LHC limits. Further, for electroweak gauginos and sleptons, the 
current LHC limits will also be weakened in our case for the same reason. After that we 
also require surviving samples to avoid Landau singularity at GUT scale and we checked 
that all of our surviving samples satisfy V/A? + K2 < 0.75 in NMSSM. We note that a large 
tan f exist in the surviving samples of the MSSM, this is because that a 125 GeV neutral 
Higgs mass is guaranteed by a large A, (which provides X,/M, close to V6) even for tan 8 
as large as 40. As for the flavor constraints, we projected our samples onto the tan 8 versus 
the charged Higgs mass plane and found that when tan f increases the charged Higgs mass 
grows dramatically (especially, for tan 8 close to 40, the charged Higgs mass is heavier than 
700 GeV) and thus can satisfy the flavor constraints. For the samples surviving the above 
constraints (i)-(iv), we further perform a fit by using the available Higgs data at the LHC. 
We define the Higgs signal strength ju; as 


i 
X ope Br; 


< NUN 18 
H 30748 B^ ( ) 


where p is the Higgs boson production mode and 7 stands for the measured channels by 
Tevatron, ATLAS and CMS collaborations. For each production mode p, its contribution 


to the channel 7 can be determined by the selection efficiency é, 5g. We summarize all 


experimental signal strength y;*? with their 1o error-bars and selection efficiencies in Fig [I] 
We can see that most measurement results are consistent with the SM predictions. The 
CMS and ATLAS collaborations also reported their observations of the Higgs mass M; 


| 
125.8 + 0.5 + 0.2 GeV (CMS ZZ), 


( 
125.4 + 0.5 + 0.6 GeV (CMS yy), 
124.3 + 0.6 + 0.5 GeV (ATLAS ZZ), 
126.8 + 0.2 + 0.7 GeV (ATLAS ^). 


Mg = (19) 


We use the combined Higgs mass M,"” = 125.66 + 0.34 cevied. The x? definition in our 
fit is 


22 
(u—u y? MMT 
i=1 O; OM, 


where c; and cy, only denote the experimental errors. 


B. The cross section of bb — hh with SUSY-QCD correction 


We use CTEQ6L1 and CTEQ6m for the leading order and SUSY-QCD calculation, 


The renormalization scale upg and factorization scale up basically can vary 


respectively. 


between M,/2 and 2M,. In order to compare our results with where ug = up = My/2 


is assumed, we also made this assumption in our calculation. The input parameters of the 


SM are taken as |62 


1000 


my = 4.7 GeV, m, = 173.1 GeV, mz 


— 91.19 GeV, 


sin? Oy = 0.2228, a,(m,) = 0.1033, a = 1/128. (21) 
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FIG. 2: The scatter plot of the parameter space satisfying the experimental constraints (i-iv), 


showing the hadronic cross sections of the SM-like Higgs pair productions via bb annihilation (with 


SUSY QCD correction) and gg fusion versus M4 and tan f at the 14 TeV LHC in MSSM and 


NMSSM. 


In Fig. 2] we display the parameter space satisfying the experimental constraints (i-iv), 
showing the cross sections of the SM-like Higgs pair productions via bb annihilation (with 
SUSY QCD correction) and gg fusion versus M4 at the 14 TeV LHC in MSSM and NMSSM. 
In this paper we aim to investigate the property of the bb — hh production by including 
the SUSY QCD corrections. For the gg — hh production, we only calculate its cross 
section at one-loop level, not including the SUSY QCD corrections due to its small relative 
correction 63 comparing the SUSY QCD correction on bb — hh process. We used our own 
codes and combined them with Looptools to do our calculation. We checked our results 
with [27] and found good agreement. 

We checked that our results agree vith[p9 for bb — hh and with for the gluon 
fusion process. We can see that due to the constraints from the LHC and B-physics, such as 
H/A > rr and B, > u*u- lcd. the values of m4 must be larger than about 300 
GeV. In the MSSM the maximal cross section can still reach 50 fb at 14 TeV LHC, which 
can be competitive with gg — hh. However, we also notice that the hadronic cross section 
proceeding through bb + hh deceases when ma or tan 8 becomes large. The reason can be 
understood as follows. On the one hand, for a moderate m4, the dominant contribution to 
bb — hh comes from the resonant production bb + H — hh. With the increase of M4, the 
mass of H gets heavy and then the production rate of bb — hh is suppressed. Besides, the 
coupling of hhH will approach to zero for a large m4 and also leads to the reduction of the 
cross section. On the other hand, for a small tan 8, H has a large branching ratio into a pair 
of Higgses h^|65], for a large tan 3, the production rate of bb > H can be enhanced but the 
branch ratio of H — hh is highly suppressed. So the total production rate of bb > hh will 
become small.The decoupling behavior of the cross section proceeding through gg — hh can 
be understood with the following considerations: To predict a 125 Gev Higgs boson, a large 
A, is required, which induces a sizable SUSY effect for the process gg — hh. Ma affects 
the process gg — hh mainly through the Higgs mass mj. So when we require m; in the 
range of 123-127 GeV, the process gg — hh is not sensitive to MA. Further, since gg > hh 
is dominated by the stop loops, the value of tan 8 affects this process through the coupling 
ht it;. Because this coupling is not sensitive to tan 8 for our surviving points, our results 
depend weakly on tan f. 

In NMSSM the SM-like Higgs boson h with mass around 125 GeV can be either hı or ho. 


However, we focus on the h = hg scenario that is more welcomed by the naturalness. From 
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FigD]we can see that the maximal cross section of bb — hh can only reach about 4 fb, which 
is much smaller than gg — hh. We find that the suppression of bb — hh in NMSSM mainly 
has two reasons. One is that in NMSSM the tan @ value is around 3-5 which is much smaller 
than in MSSM which is always larger than 10. So the tan 8 enhancement on h,bb coupling 
is not significant in NMSSM. The other reason is the h3h3h3 coupling is suppressed for most 
surviving points (the main reason is the cancelation of different contributions). Besides, in 
the NMSSM the 125 GeV Higgs mass requires a small tan and a large A. So the cross 
section of bb — hh can hardly enhanced by tan £. 
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M (GeV) tanf 


600 800 1000 1000 1500 2000 2500 3000 
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FIG. 3: Same as Fig.2, but showing the relative SUSY-QCD correction for the bb — Ah in the 
MSSM. Here the samples satisfying the experimental constraints (i-iv) are further classified ac- 


cording to the Higgs data: within lo (red dots), outside 1e but within 2c (green triangles). 


To further investigate the influence of the Higgs data in Fig.2 on the SUSY-QCD effect 
in bb — hh, we define the relative SUSY-QCD correction dsacp as 


OSQCD — OLO (22) 


ÓsQCD = 
OLO 
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In our calculation we use the oZ? for the LO cross-section and aN”? for the NLO cross- 
sections, respectively. In Fig. [3] we show the dependence of dsgcp for the bb — hh on the 
SUSY parameters M4, tan f, the lightest sbottom mass (mj, ) and gluino mass (mg) in the 
MSSM. In this figure the samples satisfying the experimental constraints (i-iv) are further 
classified according to the Higgs data: We use the x? and the degree of freedom to calculate 
the p-value for each point and plot the points whose p-values are larger than 0.045 (2c) and 
0.318 (1e). The degree of freedom is 15 [23(experimental observables)-8(free parameters)| 
for MSSM and 12 [23(experimental observables)-11(free parameters)| for NMSSM.From the 
upper panel we can see that a heavy m4 (> 400 GeV) and a moderate tan 8 (10 ~ 40) are 
favored by the Higgs data and the SUSY-QCD correction can maximally reach about 4596 
for the samples in lo range. Similar to Fig[2] dsgcp decreases when m4 becomes heavy. 
From the lower panel we note that for heavy m;, and mg, the SUSY-QCD effects decouple 
slowly. This behavior is because that the SUSY-QCD corrections depend on the ratio of 
the SUSY parameters. For example, in the triangle diagrams, the SUSY-QCD correction 
to the vertex hbb is proportional to M2,,/M2 and Maw / M? ad kd. So only when all 
the sparticles and m4 are heavy, the SUSY-QCD effect can completely decouple from the 
process of bb — hh. 

The relative SUSY-QCD corrections for the bb — Ah in the NMSSM are presented in 
Fig] It can be seen that the maximal SUSY-QCD correction can reach 15% for the samples 
in lo range. From the upper panel we can see that dsgcp becomes small with the increase 
of A or Maa. The reason is that with the increase of the A, the mp, gets heavy and its 
contribution to the cross section becomes small. From the lower panel we see that, due to 
the residual effects of the sparticles, the SUSY-QCD corrections can still reach about 9% 
for heavy sbottom and gluino. 

In Fig[5] we show the total cross section of the Higgs pair production at the 14 TeV LHC 
(via both bb annihilation and gg fusion) for the samples in the 1c and 2c ranges of the Higgs 
data. We can see that in the lo range the total cross section can be maximally enhanced 
by a factor of 2.7 and 2.2 in the MSSM and NMSSM, respectively. 

Finally, considering the null results of the direct search for sparticles at the LHC, we 
investigate the SUSY-QCD effect in Higgs pair production in the limit of heavy sparticles. 
For simplicity, we assume a common mass Msysy for all relevant SUSY mass parameters: 


Msusy = Mg = Mp = A; = Ay = M; = M,. In FigÉwe display the ratio uoc p 
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FIG. 4: Same as FigB] but for the NMSSM. 


We can see that for Msusy = 1 TeV, the ratios will maximally reach 3 and 2 in the MSSM 
and NMSSM, respectively. When Msysy goes up to 5 TeV, the enhancements become weak 
but can still reach 1.8 and 1.4 in the MSSM and NMSSM, respectively. So the effects of 
heavy sparticles decouple quite slowly from the Higgs pair production. We checked that the 
SUSY effects decouple quickly in bb — hh but slowly in gg — hh. 


IV. CONCLUSION 


We considered the current experimental constraints on the parameter space of the MSSM 
and NMSSM. Then in the allowed parameter space we examined bb — hh (h is the 125 GeV 
SM-like Higg boson) with one-loop SUSY QCD correction and compared it with gg — hh. 
We obtained the following observations: (i) For the MSSM the production rate of bb > hh 
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FIG. 5: The total cross section of the Higgs pair production at the 14 TeV LHC via both bb 
annihilation (include the SQCD correction) and gg fusion (without the SQCD correction) in MSSM 


and NMSSM. 


(with one-loop SUSY QCD correction) can reach 50 fb and thus can be competitive with 
gg — hh, while for the NMSSM bb — hh has a much smaller rate than gg — hh due to 
the suppression of the Abb coupling ; (ii) The SUSY-QCD correction to bb — hh is sizable, 
which can reach 45% for the MSSM and 15% for the NMSSM within the lo region of the 
Higgs data; (ii) In the heavy SUSY limit (all soft mass parameters become heavy), the 
SUSY effects decouple rather slowly from the Higgs pair production, which, for Msusy = 5 
TeV and my < 1 TeV, can enhance the production rate by a factor of 1.5 and 1.3 for the 
MSSM and NMSSM, respectively. Therefore, the Higgs pair production may be helpful for 
unraveling the effects of heavy SUSY. 
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FIG. 6: The cross section of Higgs pair production via both bb annihilation (include the SQCD 
correction) and gg fusion (without the SQCD correction) in MSSM and NMSSM for heavy sparticle 
masses at 14 TeV LHC. 
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